The exact solution of the eigenproblem for the 
parametric down conversion process in the Kerr 

medium 

Goce Chadzitaskos*, Maciej Horowski**, Anatol Odzijewicz**, Agnieszka 

Tereszkiewicz** ^ 



* Faculty of Nuclear sciences and Physical Engineering 
CTU, Brehova 7, CZ-115 19 Praha 1, Czech Repubhc 
**Institute of Theoretical Physics 
University in Bialystok 
Lipowa 41, 15-424 Bialystok, Poland 



Abstract 

The eigenproblem for a class of Hamiltonians of the parametric down conver- 
sion process in the Kerr medium is solved. Some physical characteristics of the 
system are calculated. 

1 Introduction 

During the recent years the parametric down conversion process has been intensively 
investigated. It is the simplest nonlinear model in which entangled photons (biphotons) 
are produced, see |P-Lj . Another intensively studied phenomenon in quantum optics 



is the Kerr effect, which is the simplest example of a optical bistability, |P-V-Lj . In 
this paper we consider the Hamiltonian which describes both effects simultaneously. 
Following llj we chose the special form of the parametric down conversion term and 
therefore our Hamiltonian is given by ()2.3|1 . Using the techniques developed in |()-H-T| 
IH-O-Tj . which are strictly related to the theory of orthogonal polynomials, we solve 
the eigenproblem for Hamiltonian ()2.3p . This method works only if the parameters 
Ui, U2, g and Ki, K2 included in the Hamiltonian, satisfy extra conditions ()3.8p . and 
then the solution is given in terms of the Dual Hahn polynomials. In the last section 
we show the way how to calculate the expecattion values of physical quantities in any 
state, in particular coherent state and occupation number state. Physical importance 
of this model required to take in consideration more details, and it will be done in a 
separate paper. 

^E-mail:**chadzita@br. fjfi.cvut.cz, *aodzijew@labfiz. uwb.edu. pi, horowski@alpha.uwb.edu.pl, 
tereszk@alpha.uwb.edu.pl 
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2 The Hamiltonian of the system 

Let a^^, ci'iid a^, be the usual annihilation and creation operators of photons in the 
mode 1 and 2. They satisfy the standard commutation relations 

[a„a*] = 6^j, [a,,aj] = 0, i,j = 1,2, (2.1) 

and act in the Hilbert space H spanned by the Fock basis 

Bf ■■= I \ni, 712) = ^ ' |0, 0) : ni,n2 = 0,l,...} (2.2) 

in the usual way. Process of quantum-optical parametric down conversion is described 
in any Hamiltonian by the term afa2, if additionally we need an intensity dependent 
process we must multiply this term by a function of a^a^ and a2a2. In this paper, due 
to we consider the system with the amplitude dependent down conversion process 
described by g^/a^a^ afag, where g is a real constant. The Kerr effect is described by 
the terms of the form (a*aj)^. Thus in order to take into consideration the free energy 
of the modes we will study the following Hamiltonian 

H = uiala, + uj2a;a2 + Ki{ala,y + K2(a*2a2)' + g a?a* + af a2 (2.3) 

where ui, U2, Ki, K2 are real constants. It is easy to check that the operator 

R := a^a^ + 2a;a2 (2.4) 

and the projective operator P defined by 

I \ / if '^i is even , . 

1 77-1, 712) it Til IS odd 

are constants of motion, and each vector of Fock basis Bp is its eigenvector. In order 
to decompose the Hilbert space Ti into invariant subspaces let us notice that for each 
nonegative integers M and p E {0,1} the subspace 

Hp,M := span{|2A; + j9,M- /c), A; = 0, 1, . . . , M} (2.6) 
is the eigenspace of the operators R and P, i.e. 

^\iP)E7ip,M Yi\ilj) = {2M + p)\iP) V\'ilj)=p\^). (2.7) 
Because of this fact the original Fock space can be decomposed 

^= ^P,M (2.8) 

PS{0,1} 
Me N U{0} 

where, of course , each Tip^M is H- invariant and of finite dimension 

dimHp,M = M + 1. (2.9) 
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The eigenproblem for the Hamiltonian ()2.3p is reduced to the eigenproblems of the 
finite-dimensional operators Hpjvf := HIt^^^j-. The matrix form of iip,M in the basis 



{\2k + p,M - 


















/«o bo 





■■■ \ 




A 








• •■ o\ 






bo ai 


bi 









1 











61 


02 













1 


■•• 


, (2.10) 








bM-1 










■•• 






[0 ■■■ 





bu-i clm / 




VO 







1) 





where 



bk 



k 



1, 



(4Ki + K2) + — (2cui -UJ2 + ApK, - 2MK2) + M{p + 
Ig Ig Z 

L 



{M-k)\l{k + l){k + p+-), 



(2.11) 



and the constant is 

Cp,Af 



-2gM{p+- 



ujip + UJ2M + /Cip + K2M^. 



(2.12) 



In the next section we show that under some special conditions for parameters uji, 102^ 
g, Ki, K2 the first matrix of ()2.10|) can be diagonalized by the Dual Hahn polynomials 
theory. 



3 Dual Hahn polynomials and spectrum of the Hamil- 
tonian 

We start this section with the brief definition and some important properties of the 
Dual Hahn polynomials |K-Sj . For any fixed natural number and for 7>— 1,5>— 1 
let us define a finite set of real points 

A; = /(/ + 7 + 5 + 1), / = 0,1,...,A^ . (3.1) 

The Dual Hahn orthonormal polynomials are, by definition, the finite family of the 
polynomials {Pk{Xu 7, 5, N)}^^^ of the discrete variable A/ satisfying the following three 
term recurrence formula 

\iPk{\u 7, 5. N) = bkPk+i{\u 7, <5, N) + akPk{\u 7, '5, N) + hk^iPk-i{\u 7, N), (3.2) 

where a^, hk are defined by 

ak = {k + j + l){N-k) + k{N + 6 + l-k) 

(3.3) 

bk = ^/{k + l){k + J + 1){N - k){N + 6 - k). 
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Formula ()3.2|) must be satisfied for any k,l = 0,1, . . . , N. The ortlionormality means 

N 

J2 Qn{1; 1, S)Pr{Xi; 7, S, N)Ps{Xf, 7, 6, N) = 5rs, (3.4) 
1=0 

wfiere 

(-!)'(/ + 7 + d + l)Ar+l(() + l)i II 

and tlie Pocliliammer-symbol {a)i is defined by tlie gamma function 

r(a + /) 

a i := . 3.6 

r(a) 

Let us note tfiat tlie Dual Halm orthonormal polynomials can be expressed in terms of 
the hypergeometric functions 



Now we use this polynomials to solve the eigenproblem of the Hamiltonian ()2.1U|) for 
the case when ui, U2, g, Ki, K2 obey special conditions 

2ijji — 002 — g = 

2Ki + g = (3.8) 
K2 + 2g = 0. 

We observe that the substitution 

7 = p-^ 5 = N = M (3.9) 

realizes the equivalence between the formulae ()2.1H1 and Hence the eigenproblem 

Hp,Af |^i,p,Af) = Ei^p^m\Ei^p^m) I = 0,1, ... M (3.10) 
due to (j3.2j) is solved in terms of Dual Hahn polynomials, 

Ei^p^M = 2gXi-2gM(^p+^^+u;ip + L02M + Kip^ + K2M^ (3.11) 

M 



Ei,p,m) = 'iL/(/;p-i 0)|^Pfe(A/;p-^,0;M) \2k + p,M-k). (3.12) 



2 

fc=0 

The multiplier on the right hand side of the last formula is chosen to fulfil the orthonor- 
mality condition 

{Ek,p,M\Ek',p',M') = 5kk' Spp' Sma4', (3.13) 

what follows from ()3.4j) . We can also express the vectors \2k + p, M — k) in the basis 
of the eigenvectors of Hp_ 



p,M 



1=0 



*^ I i — 1 

2k + p,M-k) := J2JgM{l;p-7^,0)Pk{Xr,p--,0;M) \Ei^p,m). (3.14) 
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4 Physical characteristics of the system 

In this section we describe the time evolution of expectation values of physically inter- 
esting observable on arbitrary state To do this, it is important to calculate the 
matrix elements of the evolution operator e^*^*. Because Tip^u is H invariant, the only 
nonvanishing matrix elements of e~*^* are of the form 

{2l + p,M -l\e-'^^ \2k + p,M -k) = (4.1) 

^ 111 
= ^ e-* * ^...,M ^^^(^•. p__^ 0)Pi{X^;p - -, 0, M)PkiXr,p - -, 0, M) 

what follows from ()3.14|) and ()3.10p . The decomposition ()2.8|) leads us to the following 
resolution of identity operator 1 in 7i 

oo M 1 

1=^^"^ \2k + p,M -k){2k + p,M -k\. (4.2) 

M=0 k=0 p=0 

Let X be any observable and any arbitrary chosen initial state. Then using four 
times the resolution ()4.2p and ()4.H) we obtain that the time evolution of the expectation 
value of X in the state 1-0) is given by 

oo 1 Ml M4 

(^|e'H*Xe-«*|V^) = J2 J2 H 5^ (V^|2A:i+pi,Mi-fci)x (4.3) 

Mi,M4=0 pi,P4=0 fci,fc2=0 fc3,fc4=0 

X ( 2A;i + pi. Ml - ki |e*"* \2k2 + pi. Mi - A;2 ) x 

X {2k2 +Pi,Mi-k2\X2k3+ p4, M4 - ks) X 

x( 2A;3 +p4,M4 - A;3 |e"*"*|fc4 + p4, M4 - k^ ) (A;4 + p4,M4 - k^l tp). 
Substituting = \2l + p, M — I) in ()4.3|) and using ()4.ip we obtain as a finite sum 
(2/ + p,M-/|e^"*Xe-^"*|2/ + p,M-/) = (4.4) 

M 

= e-2^*5(^^^"^^i)PKA,,;p--,0,M)P,,(A,,;p--,0,M) x 

111 1 

x^M(ii;p- -,0) gMU2;p- ^,0)Pks{Xj^;p- -,0,M)Pi{\j^;p- -,0,M) X 
x{2k2+p,M - k2\ X \2k3+p,M - ks). 

Moreover, when we put for example X(t) = ai(t)ai(t), then 

(2/+p,M-/|e*"*a*(t)ai(t)e-^^*|2/ + j9,M-/) = (4.5) 

AI 

= p + 2 J2 ke-^''^^^^^-^^^^ Pi{X,,;p--,0,M)Pk{X,,;p--,0,M)x 

ji,j2,k=0 

XQM{ji]p-^,0) gM{32]p-^,0)Pk{Xj^]p-^,0,M)Pi{Xj^;p-^,0,M). 

Now the calculation of the dispersion for this operator is straightforward, too. When 
we consider a coherent state lip) = \zi, Z2) = n2=o v'W^n2! I^i' '"'2) then ()4.3p takes 
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the form 



oo 1 Ml M4 ^Pi+P4 3fi+M4 /^2\ki+ki 

M,^=opiir=o.,tlo.37r=o V(2fci+Pi)!(2fc4+J>4)!(Mi-/c0!(M4-A:4)! V^2, 
x(2A;i+pi,Mi-fci |e^"* |2A;2 + pi, Mi - ^2 ) x 

X (2A;2 + pi, Ml - A;2| X 2A;3 + p^, M4 - A;3) ( 2h + P4, M4 - h \e-'^' \2ki + p4, M4 - 

The formulas ()4.3|) - ()4.6p initiahze a new area of investigations and numeric calculations. 
It is another thing to study this ideas according to the spirit of the theory presented 
in |H-()-Tj . This short paper do not allow us to expand this subject. Because of the 
physical importance of this theme we will devote to this subject a separate paper. 
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